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Abstract. We study the effect on the stationary currents of constraints affecting the hopping rates in stochastic 
particle systems. In the framework of Zero Range Processes with drift within a finite volume, we discuss 
how the current is reduced by the presence of the constraint and deduce exact formulae, fully explicit in 
some cases. The model discussed here has been introduced in Ref. [1] and is relevant for the description of 
pedestrian motion in elongated dark corridors, where the constraint on the hopping rates can be related to 
limitations on the interaction distance among pedestrians. 
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1. Introduction 

This paper reports on exact results for the calculation of stationary currents for a class of one-dimensional 
zero-range processes with threshold modeling the dynamics of pedestrians walking in an elongated corridor 
with no visibility. Modeling the corridor as a one-dimensional array of discrete sites, we assume that more 
pedestrians (particles) can occupy the same site (forming, possibly, social structures) and no interaction 
between these particles takes place. 

The evolution of the pedestrians is determined by the level of occupancy of the sites. The main specific 
feature is the presence of the activation threshold which keeps the “escape rate” minimal until a certain 
occupation number on the site, corresponding to the threshold, is reached. The threshold can be related to 
the limited interaction distance among pedestrians Ref. |l] : only if a site is sufficiently populated pedestrians 
can efficiently exchange information and move coherently to a neighboring spot. The approach can be further 
extended to consider the presence of multiple thresholds (e.g. communication saturation thresholds cf. Ref. 
[2] , de-centralized task-allocation thresholds cf. Ref. [3], and so on). However, in that case of multiple 
thresholds exact calculations are out of reach. Our attempt is particularly relevant for the construction of 
exact microscopic and macroscopic fundamental diagrams (explicit relationships between the pedestrians 
speed and local density, see Ref. m for pedestrians motion in one-dimensional models. 

The distinguishing feature of the model is the presence of the activation threshold whose meaning for 
pedestrian motion has been discussed above (see, also, Ref. [5] for the discussion of threshold effects in 
pedestrian dynamics in the framework of a two-dimensional model). Nevertheless, different interpretations 
are possible: for instance, in the framework of Porous Media , the bulk porosity estimates how many particles 
can be accommodated in a cell and this connects to the activation threshold. Activation thresholds are also 
meaningful in pure mechanical applications: imagine that a device is equipped with valve-like door whose 
opening results from the balance between the pressure inside the cell and an external force exerted by a 
spring. A minimal - structural - opening of the door, with the spring maintained at rest, corresponds to the 
presence of an activation threshold. A psychologico-geometrical interpretation is also possible: the activation 
threshold can be indeed regarded as a measure of the domain of communication between the individuals and 
the level this communication is processed towards a decision on the motion. 
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Moreover, the mathematical framework developed in this work paves also the way for a deeper understand¬ 
ing, through the prism of stochastic dynamics, of the kinetic mechanisms giving rise to the hydrodynamic 
properties observed, e.g., in the study of the transport of a gas or a liquid through polymeric matrices, see 
e.g. nm. Interestingly, note that despite the microscopic dynamics described in the sequel is not related 
to any energy function, the activation threshold present in our model connects, quite naturally, with the 
activation energy occurring in the Arrhenius expression for the rate of a chemical reaction, with the site 
occupancy (a random variable) playing the role of a temperature. 

It is also worth mentioning that a suitable variant of the model discussed below was also introduced in the 
literature, see Ref. [5] , to investigate the thermodynamic properties of heterogeneous materials, in which, 
e.g., a single site may be equipped with a hopping rate whose dependence on the site occupancy differs 
from the rule assigned to the remaining sites. This situation was shown to give rise to interesting physical 
phenomena, cf. also Ref. [9]. 

Coming back to the original problem, the microscopic dynamics is modeled here via a Zero Range Process 
(ZRP), cf. Ref. m, i n which the particles hop, with a certain intensity and an assigned probability, to the 
neighboring sites and in which the threshold affects the intensity of the jumps from each lattice site. In 
the framework of ZRP models, thresholds are not a novelty, see e.g. Refs. EH 23 where condensation and 
metastability effects have been studied. In those papers the value of the threshold is scaled with the size of 
the system and distinguishes between “fast” sites, namely those with a sufficiently small number of particles, 
and “slow” sites, the remaining ones. 

We exploit the threshold in a different fashion, see Ref. El 13: indeed, for our application, the hopping 
rate must be increasing with the number of particles on the spot and the threshold is used to activate the 
regime in which the rate starts to increase linearly with the number of particles. In Ref. El 13 the model 
has been studied in the hydrodynamic limit, whereas in this Note we solve the model for finite values of the 
lattice size and the number of particles. In particular, here we compute the steady state current, which is, 
even in the pedestrian motion interpretation, the main quantity of interest. 

In the absence of threshold, the stationary current increases proportionally to the number of particles, 
whereas it tends to an asymptotic value when the threshold is equal to the number of particles. In such a 
case no site exceeds the threshold and the hopping rate stays always equal to its minimal value. We compute 
the steady current for any intermediate value of the threshold and, in particular, we prove that thresholds 
proportional to the number of particles are sufficient to induce the asymptotic limiting regime. 

We remark that the focus of the paper is on the combined effect of a bias (i.e., a driving force, breaking 
the condition of detailed balance) and an activation threshold in presence of a finite number of particles 
moving on a finite lattice , endowed with periodic boundary conditions. Thus, we shed light on the finite size 
corrections to the value of the stationary current obtained in the hydrodynamic limit of the model (see Ref. 
E3 for mathematical details): this program is pursued, here, by evaluating the canonical partition function, 
which can be explicitly read out in a few cases. 

The paper is organized as follows. In Section [2] we introduce the model and define the stationary current. 
In Section [3] we derive the expression of the partition function of the model that is exploited in Section [4] to 
compute the current and to compare theoretical results to numerical simulations. 

Conclusions are drawn in Section [5] 


2. The model 

We consider a positive integer L and define a ZRP on the finite torus A := {1,..., L} C Z, cf. Refs. 1181 ' 141 . 
We fix N £ Z + and consider the finite state space Ql,n'- 

l 

n L ,N = {v G {0,..., 7V} A , £ r h: = IV } . (2.1) 

X — l 

Given rj = ( 771 ,..., tjl) £ 11 l,n the integer p x is called number of particle at site x £ A in the state or 
configuration p. Pick the threshold T £ {1,..., N} and define the intensity 

gr{ 0) = 0, gr(k) = 1 for 1 < k <T, gr{k) = k — T + 1 for k > T . (2-2) 
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The ZRP considered in this paper is the continuous time Markov process rj t £ l,n, with t > 0, such 
that each site a; £ A is updated with intensity gr(Vx(t)) and, once such a site x is chosen, a particle jumps to 
the neighboring sites x — 1 and x + 1 with probabilities, respectively, 1 — p and p (recall periodic boundary 
conditions are imposed). Note that the equilibrium condition of detailed balance holds only for p = 1/2. 
The above described jump process corresponds, hence, to an inhomogeneous Poisson process with hopping 
rates 

r x ,x+i(v) = g(Vx)p, r XiX _ 1 {if) = g(rj x )(l - p) and r XtV {rj) = 0 for y ± x - 1, x + 1 . (2.3) 


Given the threshold T, the intensity function is constantly equal to one up to T and then it increases 
linearly with the number of particles occupying the site. In other words, all sites with number of particles 
smaller or equal to T are treated equally by the dynamics, whereas the updating of those sites with more 
than T particles is favored. For this reason T is called activation threshold. 

We note that in the limiting case T = 1 the intensity function becomes g\{k) = fc, for k > 0, and thus 
the well known independent particle model is recovered. A different limiting situation is the one in which 
the intensity function is constantly equal to 1 for any k > 1 and equal to zero for k = 0. In this case a Zero 
Range process whose configurations can be mapped to the simple exclusion model states is found. We shall 
refer to the latter case as to the simple exclusion-Wke model. Such a model is found, in our set-up, when 
T = N . We stress that one of the interesting issues of our model is the fact that it is able to tune between 
two very different dynamics: the independent particle and simple exclusion- like behavior. 

It can be proven (see Ref. Duma) that the invariant measure of the ZRP process is a product measure 
of the form 




1 

Zl,n,t 


n 

x=l 


1 

gr{ i)• --gTivx) 


(2.4) 


for any p £ where the partition function Z^jyx is the normalization constant. 

The main quantity of interest, in our study, is the stationary current representing the difference between 
the average number of particles crossing a bond between two given sites from the left to the right and that 
in the opposite direction. More precisely, since periodic boundary conditions are imposed, the current does 
not depend on the chosen bond and is defined as 


Jl,n,t ■= {r x , x +i - r x +i, x ) l,n,t = ( 2 p— l)(gr)L,N,T 


(2.5) 


where we introduced the notation ( f(p))L,N,T '■= '52 irj k L L,N,T{'g)f( r D f° r an Y function / : £Il,n,t K- 

A general expression for the expectation, with respect to the invariant measure (2.4), of the intensity 
function can be provided (see Ref. my More precisely it holds 


(9t)l,n,t — 


Zl,n-i,t 

Zl,n,t 


( 2 . 6 ) 


Equations (2.5) and (2.6) yield the following expression 


Jl,n,t = (2 p — 1) 


Zl,n-i,t 

Zl,n,t 


(2.7) 


for steady state current. 


3. Canonical partition function 

The final goal of this paper is computing the steady state current at finite volume for any value of the 
threshold. In order to apply equation (2.7) we need an explicit expression of the partition function. 

In this Section we shall prove an exact formula expressing the partition function in terms of sums of 
factorials and yielding explicit expression of the partition function in the limiting cases T = 1 and T = N. 

We first state a combinatorial lemma whose proof is based on techniques borrowed from Ref. [T5] . Given 
the positive integers i,j-,k 7 we let <3?(i,j, k) be the number of ways in which j indistinguishable balls can be 
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distributed into i distinguishable urns with at most k balls into each urn. Note that for j > ki we shall 
understand k) = 0. For i,j positive integers, we also let 



(3.8) 


which can be proven to be equal to the number of ways in which j indistinguishable balls can be distributed 
into i distinguishable urns, see Ref. m section 3.2.12]. 

Lemma 3.1. Let i,j , k positive integers such that j < ki, then 


\ ■ _ 

s—o 


i(k + 1) 


(3.9) 


where s := min{z, [j/(k + 1)J}. 


We omit the simple proof of the Lemma (3.11. It suffices to remark, here, that the proof relies on the 
theory of generating functions, as presented, e.g., in Ref. [TS] section 3.3.2]. 

The expression of <I>(i,j, k) provided by the Lemma (3.1) attains a simpler form in a few cases. For 
instance, when j = k , no constraint is imposed on the allocation of balls among the urns, hence one should 
find 

= (3-10) 

This is indeed the case, since it holds s = 0. Note that this is the result which is found when, in the 
Bose-Einstein statistics, one counts the number of ways in which j particles can be distributed among i 
states. A second relevant case is the one in which at most one particle can be allocated into each urn. The 
corresponding value of $ may then be derived either from Eq. (3.9), by using the fact that, since k = 1 and 
i > j , it holds s = [j/2 \, or from the combinatorial definition of <f>(i, j, 1). In either case, one obtains 


®(hj, 1) = 


(3.11) 


Note that this is the result one encounters in the Fermi-Dirac statistics, in which one counts the number 
of ways in which j particles can be distributed among i states with the limitation, due to the exclusion 
principle, of at most one particle per state. 

We can now state our main result about the canonical partition function of the ZRP model. Recall, see 
Eq. (|2.4|), that 

1 (3-12) 


z l ,n,t = yi n 


77: 1 77I =iV x —1 


g T (l) ■ ■ ■ g T (Vx) 


Theorem 3.2. For N, L positive integers 


z - LN 


Moreover, for any T > 2 and L > \N/(T — 1)] 


(3.13) 


Z l , n , t = HL,N,T-2)+Y / 


m 


J2®(L-m,N- [(T- l)m + n],T - 2)^- (3.14) 

n —0 


where m = \N/{T — 1)J and n(m) = N — (T — 1 )m. 

Proof. Consider, first, the case T = 1. By (2.21 and ( |3.12[ ) we get 

z l ,n,i= e n “i 


77:177I = AT x=l 
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where we also used the convention 0! = 1. Thus, equation (3.131 follows immediately by applying the 
multinomial theorem, see, e.g., Ref. B33 equation (3.35)]. 

Consider, now, the case T > 2. Call m £ [0, to] the number of sites in which the number of particles is 
larger than T — 1 and n £ [0, n] the number of particles that, for a given to, exceeds the value T. Given a 
configuration rj, let also 


n x = 0 for i) x <T — 1 and n x = rj x — (T — 1) for r) x > T ~ 1 . 
Then, the partition function can be rewritten as 


(3.15) 


m 


n 


Z LjN: t = HL,N,T- 2)+ ^ [\}_^<f>(L-m,N-[(T-l)m + n],T-2) 

m—1 ' ' n—0 


z^ 

n x >0 


ni! ... n m ! 


The first term in Eq. ^ takes into account the contribution to the sum defining the partition function of 
those configurations in which no site has a number of particles larger or equal to T — 1. The second term 
can be explained as follows: the first binomial coefficient counts the number of ways one can choose the to 
sites such that r) x > T — 1. Note that to denotes the maximum value attained by to, for which it holds: 
to = min(L, \_N/(T — 1)J). Yet, by requiring L > \N/(T — 1)], one has m = [N/(T — 1)J, as indicated in 
the statement of the Theorem. The coefficient 4>(L — m,N — \(T — 1)to + ti\,T — 2) counts the number of 
ways to allocate the remaining N — \{T — l)m + n] particles on the L — to sites for which it holds r] x < T — 2. 
The last sum counts the number of ways in which the particles exceeding T, namely, those on the top of 
the T — 1 filled columns, can be distributed on the to sites. Finally, recalling (|2.2|), we have that the last 


factor in the equation is a smart rewriting of the last factor in (3.12). Equation ( |3. 14 ) in the theorem finally 
follows by using the multinomial theorem (see, e.g., Ref. fl5l equation (3.35)]). O. 


We remark that, although (3.14) is not an explicit expression for the partition function, it is nevertheless 


very useful. Indeed, the sum over the configuration space present in the definition of the partition function, 


Eq. (3.12), involves a number of terms exponentially large in the number of particles N, whereas the sum in 


(3.14) is only polynomial in N. Moreover, the expression for partition function given in Eq. (3.12) involves 


a constraint, namely |?y| = N, which has been removed in (3.14). 

It is also interesting to remark that in the simple exclusion like regime, namely, T = N, the partition 
function Zl^n,n can be written explicitly as 


Zl,n,n = 4 >(T, N, N) = ^ 


L + N -1 
N 


(3.16) 


To prove this formula, we compute, first, the term 4 )(L,N,N — 2) in (3.14). By using (3.9) with i = L, 
j = N, and k = N — 2, noted that s = 1, one finds 


$(L, N, N - 2) = 4>(L, N, N ) - Q = 4>(L, N, N ) - L 2 


(3.17) 


Next, we evaluate the sums in Eq. ( |3. 14 ). Since to = n = 1, one needs to calculate just the terms 
4>(L — 1, 0, N — 2) and 4>(L — 1,1, N — 2). Since, for both terms, s = 0 for both terms in (3.91, we get 


4-(L - 1,0, IV- 2) = ( Q 1 )= 1 and $(£ - 1,1, N - 2) = / ^^L-l (3.18) 


Equation (3.16) finally follows from (3.17), (3.18), and (3.14). 

Note that the result in Eq. (3.16) could also be directly deduced by equation fl3T2l ). Indeed, from the 
definition (2.2) of the intensity function, it follows that, for T = N, the sum in equation (3.12) is indeed a 


sum of l’s and, thus, yields straightforwardly the total number of configurations 4>(L,7V, N). 

As discussed in the Introduction, the threshold limits the hopping rate on sites whose occupancy number 
is smaller than the threshold itself, whereas, when the prescribed value of the threshold is reached, the 
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hopping rate starts increasing proportionally to the number of particles on the site. In this respect, the case 
T = TV is peculiar, because all the sites are updated with the same minimal rate regardless their occupancy 
number. 

It is also possible to guess another remarkable result: namely, when the threshold, although smaller than 
TV, scales proportionally to TV, then the stationary current is close, for large TV, to the value obtained for 
T = TV. More precisely, take a < 1, a € R + and sufficiently close to 1, and compare the canonical partition 
function of the systems with TV particles and thresholds equal, respectively, to aN and TV. 

We thus conjecture that for TV —► oo 


Zl,N,oN 

Zl,n,n 


1 + 0 ( 1 ) 


where o(l) denotes a function tending to zero in the limit TV —> oo. 

We omit, here, the lengthy algebraic details, and we just mention that this observation may be relevant 
in the study of the hydrodynamic limit of heterogeneous ZRP, in which the hopping rate from a given site 
can be modified so as to scale with the size of the system. 


4. Stationary currents and numerical simulations 

In this Section we report and compare both analytical and numerical results for the steady current in 
the ZRP with threshold introduced in Sec. [2j 

Numerics have been performed via Monte Carlo techniques by simulating the model as follows: call 
i](t) the configuration at time t, a number r is chosen at random with exponential distribution of parameter 

9T(Vx{t)) and time is update to f+r, a site is chosen at random with probability griVxit))/ XOi—i St( ? 7x(T)) 
and a particle is moved from such a site to its right with probability p and to its left with probability 1 — p. 

The Monte Carlo simulation is let, first, evolve for a number of time steps no ~ 10 7 , and the stationary 
current is thus defined as the ratio of the difference between the total number of particles jumping from site 
L to site 1 and that of particles jumping from site 1 to site L, to the total time. We remark that the initial 
number of time steps no is chosen large enough to guarantee that a constant value, with respect to time, is 
reached by the current. 

As for the analytical results on the current, note that the theory developed in the Sections above paves 
the way to the computation of the stationary current for any finite value of the parameters of the model, 

TV, L, and T. We stress that we are considering a transport problem in which a net convective flux occurs 
in the case p ^ 1/2. Equations (2.7) and (3.14) can be used to reduce the computation of the stationary 
current to an algebraic sum. In particular, in the two limiting cases T = 1 and T = TV analytic formulae can 
be derived. 

Indeed, from Eqs. (2.7) and (3.13), the steady state current for T = 1, i.e. in the independent particle 
case, reads 

TV 

Jl,n, i = (2p-l)- (4.19) 


On the other hand, Eqs. (2.7) and (3.16) imply that, for T = TV, i.e. in the simple exclusion-like regime, 
the current is given by 


Jl,n,n = (2 p - 1) 


TV 


1 




(4.20) 


We stress that the two results above are valid for any finite volume A and for any finite number of particles. 
If the limit TV, L —> oo with N/L = g is considered, the well-known hydrodynamic limit is found for the 
current, see Ref. PS equation (1.3)]. 

Coherently with its physical interpretation, the effect of the activation threshold is that of slowing down 
the current. As T is increased the steady state current decreases. In particular it is worth mentioning that 
at T = 1 the current is a linear function of the number of particles TV, whereas at T = TV the current 
saturates to a limiting value when TV is increased. For intermediate thresholds, namely, 1 < T < TV, the 
current increases slowly with TV and only after a certain value it starts growing linearly. This effect is clearly 
illustrated in Fig. 4.1 (left panel) where the current is plotted versus the total number of particle for the 
values T = 1,5,10 of the activation threshold and L = 100. 
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Figure 4.1: Current versus number of particles for L = 100. Left panel: Open circles, solid circles, and solid 
squares denote, respectively, the simulated stationary current for the threshold T = 1,5,10. The associated 
solid lines represent the analytic solutions: Eqs. (2.7) and (3.14| have been used in the cases T = 5 and 
T = 10, whereas the explicit formula in Eq. (4.19) has been used for T = 1. Right panel: Open circles and 
solid squares denote, respectively, the simulated stationary current for the threshold T = N and T = N/2. 
The associated solid lines represent the analytic solutions computed by using Eqs. (2.7) and (3.14) for 
T = N/2 and (4.20) for T = N. The two solid lines are not distinguishable in the picture. 


Data in Fig. 4.1 (right panel) refer to the cases T = N and T = N/2 , with L = 100. The saturation 
effect on the current due to the presence of the threshold T = N (simple exclusion-like regime) is clearly 
illustrated. In other words, when the hopping rate is constantly equal to one and does not depend on the 
number of particles on the site, the current tends to saturate to a constant value for N large. It is worth 
remarking that the same effect is also observed when the threshold is equal to N/ 2, suggesting that for an 
activation threshold increasing proportionally to the number of particles, the current is reduced in the same 
fashion. This property is indeed an immediate consequence of Eqs. (2.7) and (3.19). Deviation from the 
T = N behavior in the case T = N/2 can in fact be observed for small values of L and N . 


5. Conclusions 

We considered the problem of computing the steady state current in a Zero Range Process subjected to a 
drift as well as to an “activation” threshold affecting the hopping rates of the particles to the neighboring 
sites. By exploiting combinatorial arguments, we derived an exact formula for the partition function, which 
is amenable to an analytical treatment for T = 1 and T = N. We also discussed the asymptotic behavior of 
the partition function when the threshold scales proportionally to the number of particles: the latter case is 
of particular relevance in the discussion of the hydrodynamic limit of the model. We then obtained explicit 
formulae for the particle current, also supported by Monte Carlo simulations, revealing that the main effect 
of the activation threshold on the steady state dynamics is to decrease the current, thus tuning between two 
limiting regimes, the independent particle model and the simple exclusion-like process. We also remarked 
that this last behavior is shown by the model even for T < N, provided the threshold increases proportionally 
to the number of particles. 
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